Detection of Anisotropies in the Gravitational- Wave Stochastic Background 



Bruce Allen 
Department of Physics 
University of Wisconsin - Milwaukee 
PO Box 413 
Milwaukee, WI 53211, USA 
email: ballen@dirac.phys.uwm.edu 

Adrian C. Ottewill 

Department of Mathematical Physics, University College Dublin, 
Belfield, Dublin 4, Ireland 
email: ottewill@relativity.ucd.ie 
(February 7, 2008) 

By correlating the signals from a pair of gravitational-wave detectors, one can undertake sensi- 
tive searches for a stochastic background of gravitational radiation. If the stochastic background is 
anisotropic, then this correlated signal varies harmonically with the earth's rotation. We calculate 
how the harmonics of this varying signal are related to the multipole moments which characterize 
the anisotropy, and give a formula for the signal-to-noise ratio of a given harmonic. The specific 
case of the two LIGO (Laser Interferometric Gravitational Observatory) detectors, which will be- 
gin operation around the year 2000, is analyzed in detail. We consider two possible examples of 
anisotropy. If the gravitational-wave stochastic background contains a dipole intensity anisotropy 
whose origin (like that of the Cosmic Background Radiation) is motion of our local system, then 
that anisotropy will be observable by the advanced LIGO detector (with 90% confidence in one year 
of observation) if f2gw > 5.3 x 10~*/i^gg. We also study the signal produced by stochastic sources 
distributed in the same way as the luminous matter in the galactic disk, and in the same way as the 
galactic halo. The anisotropy due to sources distributed as the galactic disk or as the galactic halo 
will be observable by the advanced LIGO detector (with 90% confidence in one year of observation) 
if figw > 1.8 X 10"^°/i^oo or fig™ > 6.7 x lO'^h'^^, respectively. 



I. INTRODUCTION 



The design and construction of a number of new and more sensitive detectors of gravitational radiation is currently 
underway. These include the LIGO detector being built in the United States by a loint Caltech/MIT collaboration 
0, the VIRGO detector being built near Pisa by an Italian/French collaboration ||], the GEO-600 detector being 
built in Hannover by an Anglo/German collaboration and the TAMA-300 detector being built near Tokyo [Q. 
There are also several resonant bar detectors currently in operation, and several more refined bar and interferometric 
detectors presently in the planning and proposal stages. 

When two or more of these detectors are operating, it will become possible to correlate their signals, and in this 
way, to search for a stochastic background of gravitational radiation. The technique for such a search was originally 
described in work by Michelson ^ , Cristensen Q and Flanagan . A review of these techniques may be found in [|j . 
Such radiation might be the result of processes that took place during the very early universe. It might also result 
from the incoherent superposition of many faint unresolvable present-day sources such as coalescing binary systems. 

The stochastic gravitational-wave background might be isotropic on the sky, or it might be anisotropic. For example, 
if the background results from early-universe processes, then it might be isotropic to about the same degree as the 
2.7° K electromagnetic background radiation. On the other hand, if the background is due to white-dwarf binaries in 
our own galaxy, then they might be distributed in a pancake or bar which mimics the shape of the observed luminous 
matter in our galaxy. In this paper, we show how the correlated signal from a pair of gravitational wave detectors 
is related to multipole moments which characterize the anisotropy. This should permit a signal to be analyzed to 
search for (or place upper limits on) the multipole moments which characterize the anisotropy. In this paper, we will 
assume that the reader is already familiar with the work previously cited (references 1^-^) on stochastic background 
detection. 

These paper is organized as follows. In Section || we show how a background of stochastic gravitational radiation 
may be decomposed in a plane-wave expansion, with the coefficients of the expansion treated as stochastic random 
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variables. In Section III the properties of these random variables are related to the (frequency) spectrum and spatial 
distribution of the radiation, and a set of multipole moments are introduced which characterize the anisotropics of 
the stochastic background. These anisotropics may be searched for by studying the variations of the detector outputs 
as the earth rotates relative to the fixed cosmic frame. In Section |rV| we show how the correlation between a pair of 
detectors fixed on the earth varies with time as the earth rotates, and detail how that correlation is related to the 
anisotropics of the stochastic gravitational background. The variation of the correlation with the earth's rotation may 
be decomposed into harmonics of the earth's period. In section Section^ we introduce a set of functions 7fm(/) which 
are generalizations of the well-known overlap reduction function j{f) of references These functions characterize 

the effect of the £, m anisotropy multipole on the m'th harmonic of the detector correlation. The principal result of 
this paper is to compute these functions for the LIGO pair of detectors. This is done by introducing two special frames 



of reference, fixed with respect to the earth, in Section VI, and then performing a set of integrations in Section VII 



The final integrations are performed and explicit formula for the 7^m(/) are obtained in Section VIII. In Section IX 
we analyze the signal-to-noise ratios associated with the m'th harmonic, and give a formula which may be used to 
determine if a given anisotropy is detectable or not. Following this, we consider two specific examples of anisotropy. 
In Section ^ we consider a dipole anisotropy in the stochastic background resulting from our local proper motion. In 
order to predict the harmonics which result, it is necessary t o ad opt conventions for the normalization of the optimal 



filters; these choices are detailed in Section pCI[ In Section |XI]| we then consider the anisotropics in the stochastic 
gravitational wave background that would arise from sources distributed in the same way as our galaxy, and it's halo. 
This is followed by a short conclusion. 

Throughout this paper, c denotes the speed of light and G denotes Newton's gravitational constant. 

II. THE STOCHASTIC BACKGROUND 

The gravitational wave background may be described in terms of a perturbation to the Minkowski metric of space- 
time: 

ds^ = -c^dt^ + dx^ + habit, x)dx'^dx^. (2.1) 
In transverse traceless gauge, this can be written in the form of a plane wave expansion as 

Kb{t, x)^J2 r f '^^(/' f^)e'"^(*-^-^~/^^e^,(i7). (2.2) 

Here /ia(/, ^) is an arbitrary complex function satisfying the relation hA{—f, ^) — ^^(/, The polarization states 
are labeled hy A = +, x and is a unit vector on the two-sphere. The wave-vector of the corresponding component 
of the perturbation is fc = 2Tr ftl/c. The polarization tensors e^^ appearing in these relations may be given explicitly. 
In standard angular coordinates (0, (f) on the two-sphere one may write 

Cl = cos sin 9x + sin cf) sin 6y + cos 6z (2.3) 

m = sin 4>x — cos (py (2.4) 

n — cos cos Ox + sin cos 9y — sin 9z (2.5) 

and then choose 

6ofc(^) = ^amb - UaUb (2.6) 

^abi^) = "m-anb + nairib (2.7) 

To simplify matters later (but without any loss of generality) we assume that the z vector points along the direction 
of the earth's rotation axis. One can verify by inspection that rh and fi are a pair of orthogonal unit-length vectors 
in the plane perpendicular to O. It is simple to show that any rotation of the vectors rh and fi within the plane that 
they define simply corresponds to a trivial re-definition of the complex wave amplitudes /i+ and hx ■ 

To describe a stochastic source, we treat the complex amplitude /i^(/, f2) as a random variable with zero mean 
value. In this paper, we consider stochastic sources which are not isotropic. In principle, such a source has spectral 
properties which depends upon amplitude and frequency in an arbitrary way. For simplicity, in this paper we consider 
only stochastic sources whose directional dependence is frequency-independent. The dependence of the stochastic 
background on frequency and direction may be stated in terms of the expectation value of the product of two random 
variables /ia(/, H): 
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{h\if, n)hA'{f', n')) = SAA'Sif - ns^in, n')H{f)P{n). 



(2.8) 



Here S'^(Cl,d') is a covariant two-dimensional delta-function on the unit two-sphere. For a general stochastic source, 
the quantity H{f)P{Cl) which appears on the right hand side would be an arbitrary function of frequency and direction. 
However our assumption that the directional dependence is frequency-independent implies that the r.h.s. factors as 
shown. (Note that the expressions which we later derive for Signal-to-Noise ratios and expected signal strengths may 
be trivially extended to include the most general case.) 



III. SPECTRUM OF THE STOCHASTIC BACKGROUND 

The function H{f) determines the spectrum of the gravitational radiation. The energy density in gravitational 
waves is given by 

Pgw = ^(Aa6/l'^''), (3.1) 

where the overdot denotes a time deri vati ve, and both tensors are eval uat ed at the same space-time point {t,x). 



Substituting the plane wave expansion (2_^) into this formula and using ( p.^ ) yields 



{habit, x)h'"' {t, x)) =y] r df f dCl47r^fH{f)P{n)e^,{n)ef{n). (3.2) 
Since J2a ^ab^A ~ ^ one has 

dnP{Cl) dffH{f) = 32TT^ dnP{Cl) dffHif). (3.3) 

J-oo J JO 

In describing gravitational wave stochastic backgrounds, it is conventional to compare the energy density to the critical 
energy density Pcriticai required (today) to close the universe. This critical energy density is determined by the rate 
at which the universe is expanding today. Let us denote the Hubble expansion rate today by 

km SGC ^ 

Ho = 100 hwo —T7 = 3.2 X 10"^**/iioo sec"^ = 1.1 x lO^^^c hwo cm'^ (3.4) 

Mpc 

The value of Hq is determined by the dimensionless factor of /iioo which probably lies within the range 1/2 < /iioo < 1- 
The critical energy-density required to just close the universe is 

Pcriticai = ~ 1.6 X IQ-^hloo crgs/cm^. (3.5) 

The spectrum of an isotropic stochastic gravitational wave background is defined by a dimensionless function of 
frequency / 

!),.(/) (3.6) 
Pcriticai am/ 

Here dpgw is the energy-density in gravitational waves contained within the frequency interval (/, / + df). Using the 
definition Qgyj one obtains the relationship between the spectrum Slg^ and H{f). For / > one has 

f^gw(/) = — = f^^S2.^fHif) f dClPiCl) = ^fHif) f dnPiCl). (3.7) 
Pcriticai df Sc^H^ 32ttG J J 

This formula shows the precise interpretation of P{Cl). The stochastic background energy density is made of con- 
tributions arriving from all directions fl on the sky. The actual value of figw(/) is determined by the average value 
of P{fl); the direction dependence of this function is the same as the direction dependence of the arriving radiation 
intensity. 

For this reason, we define the multipole moments pim of the stochastic background radiation by the expansion of 
P(f2) in terms of spherical harmonic functions: 
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P(r!) = ^p,™F,™(Sl) (3.8) 

em 

where the sum is defined by 



oo oo 



E-EE - E E- (3-9) 

em e=0 m=—e m= — coe=\m\ 

In addition, without loss of generality we adopt the convention that the monopole moment is normalized by the 
condition 



Poo 



V^=> J dnP{n) = 4:TT, (3.10) 



where we assume that the spherical harmonic functions are normalized in the conventional way, so that the integrals 
of their squares over the unit sphere gives unity. Hence the spectrum of radiation is determined entirely by H{f) 
since for / > one has 



^^gw(/) = ^fH{f). (3.11) 

The directionality of the arriving radiation is determined entirely by the function P{Cl). Our fundamental assumption 
here is that the pattern of the intensity of the stochastic background is fixed in a frame of reference at rest with respect 



to the cosmological fluid. In other words, formula (3.8) for P{il) is expressed in a set of coordinates x.y,z which are 
fixed with respect to the distant stars. In those coordinates, the multipole moments pgm are constants, independent 
of time. The problem we address in this paper is this: how do we determine, from the data stream of a pair of 
interferometric detectors which are rotating with the earth, the values of (or bounds on) the multipole moments p£m? 



IV. DETECTION STRATEGY 



To determine the multipole moments pim the basic idea is to correlate the outputs of two gravitational wave 
detectors, and to look for variations of this correlated signal that are harmonics of the earth's rotational frequency. 
For this purpose, we need to consider the relationship between two different time (or frequency) scales that occur. 

The first time scale is that defined by the light travel time AT between the two sites. For the remainder of 
this section, we will assume that the two sites are the Hanford and Livingston LIGO detectors, so that AT = 
10.00 msec. The second time scale is the period of the earth's rotation about its axis relative to the cosmic frame, 
Te = 8.6 X 10* sec = 1 sidereal day. Because of the enormous disparity between these two time scales, we can define 
a third time scale, which we will refer to as the averaging time scale, r. We choose r in the range 

Ar«r<<Tc, (4.1) 

for example r = 30 sec. It is then possible to examine correlations between the two detectors as a function of time, 
averaged over periods of length r. Because t is much shorter than Tc, the correlation between the two detectors will 
vary as the earth rotates relative to the fixed cosmological frame, because of the anisotropy in P{Ct). On the other 
hand, because r is much longer than the light travel time between the two detectors, and because the detectors are 
sensitive to frequencies / « l/AT, there is a significant correlated signal on time scales shorter than r. 
Denote the output of the first detector by 

si(t) = /ii(i) + ni(t), (4.2) 

where hi is the strain due to the stochastic background and ni is the intrinsic noise of the first detector. In similar 
fashion, the output of the second detector is 

S2{t) = h2it) + n2{t). (4.3) 

Let us use the subscript i = 1, 2 to label the detectors, so for example i = 1 denotes the Hanford, WA LIGO detector 
and i = 2 denotes the Livingston, LA LIGO detector. The response hi of detector i to the gravitational radiation is 
given by 
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h,(t) = df{t)hab{t,X,{t)), 



(4.4) 



where the position of detector i's corner station is denoted by Xi{t). In this expression, the symmetric traceless tensors 
df^{t) are given by 

dfit) = I {X^m^it) - Ytm\t)) (4.5) 

where the directions of detector i's arms are defined by the unit spatial vectors X°'{t) and Yf-{t). Note that both 
d'f'{t) and Xi{t) are functions of time, because the earth rotates with respect to the cosmological rest frame. 

Define quantities which are the Fourier transforms of the signals, evaluated over an interval of one averaging time 
r centered at time t: 



j-t + T/2 

if, t) = / dt' e-^^'f'' s,{t') for z = 1, 2. (4.6) 

Jt-T/2 



These Four ier transforms are easily evaluated. Su bstituting the plane wave expansion (2.2) into the formula for the 
strain (4^) and taking the Fourier transform ( [4.6| ) we obtain 

S^{f, t)=J2 fdn r d/'e2"(/'-/)*^,(/ - f')hAif', f^)df (Oef,(f2)e-2-/'^-~'(*)/- + noise term. (4.7) 

In this expression, we have made use of the fact that the averaging time r is much less than the rotation period of 
the earth ATe, so that the vectors X°;{t), Y°'{t), and Xi{t) may be treated as constants and taken outside of the time 
integration in the Fourier transform ( |4.6| ). Wc have also defined the "finite time" approximation to the Dirac delta 
function 

^.(/)^ r'^t'e--/*'^^^, (4.8) 

J-t/2 T^f 

which reduces to the Dirac delt a fu nction 5{f) in the limit t ^ oo, but has the property that (Sr(0) = t. The final 
term on the right hand side of (4.7) is linearly proportional to the noise in detector i. 
We now define the "signal" 

/oc 
dfS\{f,t)h{f.t)Q{f) (4.9) 

where Q{J) is an optimal filter function, to be determined. Let us now determine the expectation value of S{t) 
and show how it incorporates information about the multipole moments of the stochastic background. To find the 
expected value S{t) we begin by assuming that the noise in each detector has zero mean value, and is uncorrelated 
with noise and gravitational strain in the other detector. Under these assumptions, we find 



^/Q(/)EE / '^Z' / df" / dCl / df2'e-2"W'- 

-OO A A, J — OQ J —OO J J 



•/)tg27ri(/"-/)t 



A A' 

Srif ~ f')SAf - /")<(tX(i)ea6(^^)e^j(f^')e'"^(^'^-"^^^*^-^""'-"^^^*»/^ x 

{h\{f,n)hA'ir,n')). 

We now substitute in the expectation value for the product of the amplitudes ( |2.8| ) . The integration over /" is now 
trivial. In the resulting expression, because l/r is much smaller than the "bandwidth" 1/AT of the signals, one of 
the finite-width delta functions St may be replaced by a Dirac delta function. The integration over /' is then trivial. 
The other finite- width delta function is then evaluated at zero argument, giving rise to a factor of r. One thus obtains 

{S{t)) = rdt{t)df{t) r dfQ{f)H{f) [ diiF(f2)e2-^^-^-~«/^ ^ e^,{Cl)e^,{n) (4.10) 
where Ax{t) — Xi{t) — X2{t) is the time-dependent separation vector between the two interferometer sites. 
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Not surprisingly, this previous expression can be easily simplified for the isotropic case P{(i) = 1. In this instance, 
the sum over polarizations and integral over directions can be performed explicitly, yielding (Stt/S times) a time- 
independent function of frequency known as the overlap reduction function -f{f). This overlap reduction function is 
given by 

7(/) ^ ^ dfdf J^^ dn e2"/^-A-7c (e+ (17)e+ (f2) + 6^^^^,^) • (4-11) 



Notice that in ( 4.11 ) the dependence of the positions and orientations of the detectors upon time t is not shown; this is 
because 7(/) depends only upon the relative positions and orientations, which is time (or earth-position) independent. 
Thus, in the case of an isotropic stochastic background, one finds 

P{n) = 1 => {S{t)) = —T / dfQ{f)H{f)j{f). (4.12) 

This is equation (30) of reference In the present paper, we are most interested in the anisotropic case where P{iY) 
varies with direction. In this case, the time variation of the tensors d^''{t) and Ax{t) will provide a time-dependent 
variation of the signal S{t). 

V. ROTATION HARMONICS 

Bec ause the rotation of the earth is periodic with period Tc and angular frequency We = 27r/To the expected signal 



(4.10) varies with the same period. It can therefore be represented by the Fourier series 

oo 

{S{t))= E =*. (5.1) 

m — — C30 

Because the signal is real, the amplitudes of the different harmonics satisfy Sm = S'^.^^. The amplitudes are quantities 
which would be determined by Fourier transforming the actual data: 

S,„ = ^ r dte-'"''^^'S{t). (5.2) 



Here T represents the total observation time, which later enters into ( |9.4D and which is assumed to be is a multiple of 
the earth's rotation period To. The harmonic amplitudes Sm are the (at least in principle) observable quantities on 
which any data analysis must be based; it is their expected values (Sm) which arise in the formula for the expected 
signal. Note that in an actual observation or measurement, instrumental noise in the gravitational- wave detectors 



would prevent Sm from taking on its expected value. In Section [X we analyze the typical deviations of Sm from (Sm) 
and in this way determine how accurately Sm may in fact be measured. 

Because we have assumed that the z-axis of our (cosmic) coordinate system points along the direction of the earth's 
axis, the rn'th rotation harmonic can only result from anisotropics whose phase varies with angle (j) as exp(jm0). 
These are the anisotropics associated with the Yim- Hence 

O POD °0 

(Sm) = / dfQ{f)Hif) J2 Pimllmif) (5.3) 

The functions 7£m(/) are generalizations of the overlap reduction function 7(/), which express the (frequency- 
dependent) contribution of the i?'th multipole moment to the m'th harmonic of the signal, with respect to the Earth's 
rotation. These are given by 

7^™(/) = fj dae-^"^-dt{a)df{a) j rff7 y,„(f7)e2-/^-^--(")/^ J] e;f,(f7)e^,(f2). (5.4) 

In this expression, the angle of rotation of the earth about its axis (measured from some arbitrary fiducial point) 
is denoted by a G [0, 27r) so a = LOet + constant (mod 27r). The "time-dependent" quantities d1^ and Aa; may 
equivalently be expressed as functions of a. 
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The problem at hand is now a mathematical one - to calculate the functions jemij) which are generalizations of 
the overlap reduction function 7(/). For the monopole moment = to = 0) it is easy to see that the integrand 
above is independent of earth-position a because the overlap reduction function ( |4.11 ) only depends upon the relative 
orientations of the detectors, which is a-independent, giving 



7oo(/) - (4^)-^/'7(/)- 



(5.5) 



In the next parts of this paper, we will show how to evalua te t he other 7^™- 

Our first task is to evaluate the integrals that appear in (5.4). The product df'{a)d2'^{a) is a quartic polynomial in 
sin a and cos a. One approach would be to attempt to perform the integral over O, to obtain the resulting function 
of a, and then to evaluate the integral over a. However this approach is rather cumbersome. 

A more promising method is to c ons ider t he p rojector onto the plane perpendicular to n, which may be calculated 
in terms of the vectors defined by (2^) and (2.5): 



}ab = Sab - ^a^b = TOaTOfc + naflb. 



A couple minutes of algebra starting with ( |2.6|) and (2/7) quickly establishes the identity 

^abi^)^cdi^) — QacQbd + QadQbc ~ QabQcd- 



(5.6) 



(5.7) 



We then define the set of integrals 

Cabcdia) 



dfiy,™(i))e2"*^^-^^~(")/^f)af^&Md. 



(5.8) 



The desired integrals can then be expressed in terms of this quantity. For convenience, we introduce a symbol to 
handle the contractions that occur. This is a constant tensor defined by 



(5.9) 



Making use of the fact that each of the df is symmetric in its tensor indices, and traceless, we may then write 

r2TT 



5 1 f 

limif) = / da e— <(a)d§'*(a)e^: 

STT ZTT Jq 



bcd^PgrsK^^)- 



(5.10) 



From this definition it is easy to show that 'fi^^m = {^^Y'^^li m- This follows from the parity tranformation property 

of the spherical harmonics Yirn{—^) = {^^YYem{^)- In order to now evaluate Cabcd it is convenient to introduce 
some additional coordinate systems. 



VI. COORDINATE FRAMES 



The vectors being used in this calculation are three-dimensional spatial vectors in flat Cartesian R^. Up to this 
point, we have been using a coordinate system which is fixed with respect to the cosmological fluid, and in which the 
spatial pattern of the perturbations of the stochastic background is assumed to be time-independent. This frame of 
reference is the "unprimed" frame; vectors expressed with respect to these cosmic coordinates have unprimed indices. 
We have also assumed (without any loss of generality) that the z-axis of this cosmic frame points along the direction 
of the earth's rotation axis. 

At this point, for calculational purposes, it is convenient to consider two additional coordinate systems. Thus, a 
given vector V may be expressed in terms of its components in three different frames: 

Cosmic Frame : 
Earth Frame : 
Computational Frame : V"' 

The "earth frame" is a coordinate system fixed to the earth, in which the third (z-coordinate) points along the axis 
of the earth's rotation, in the direction of the North pole. Components of vectors in this frame are denoted with 
"barred" indices. The second of these new coordinate systems will be referred to as the "calculational" coordinate 
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system. In this frame, the components of vectors are "primed". This frame is fixed with respect to the earth, and 
has its third (z-coordinate) pointing along the hne between the two gravitational-wave detectors. 

The relationship between components of vectors in these three coordinate frames may be written as matrix equations. 
Each of the matrices which appears is a special case of a rotation matrix which may be parametrized by Eulcr angles. 
Throughout this paper, we use the Euler angle conventions given by equations (4.83-6) of Afkin which arc also the 
conventions used in equations (4.5) and (4.43) of Rose jl^. It is convenient to define a pair of rotations about the z 
and y axes respectively, by 

(cos a sin a \ / cos /? — sin /3 \ 

-sina cosa Ry(/3) = 10 (6.1) 

1/ \ sin/3 cos^ / 

The most general possible rotation may be parametrized by Euler angles and is defined by the matrix R(q;, /?, 7) = 
Rz(7)Rj,(/3)R2(q;). Note that the boldface symbols here denote 3x3 square matrices. 

The matrix which relates components of vectors in the cosmic and earth frames is simply rotation through angle a 
about the z-axis: 



= i?\ = R(a, 0, 0) = R2(a) 



(6.2) 



where the first index on K labels rows and the second index labels columns, so that the operation appearing in the 
previous equation is ordinary multiplication of a column vector on the right by a square matrix on the left. Note that 
the angle a = ujet varies with time. 

Without loss of generality, assume that the freedom to choose the x- and y-axis in the Earth frame has been used 
to ensure that in this frame the separation vector Ax°' between the two detector sites has no y-component. Using 
the two LIGO sites as an example, the Earth-frame i-axis would point out from the center of the earth at an angle 
38.6881° East of the 0° line of longitude (Greenwich, England). In this frame, the coordinates of the two detector 
sites and the detector arms directions are 



Hanford, Washington : a;° = 



Livingston, Louisiana : X2 = 



707.41 ^ 
-4329.11 

4614.74 , 

3371.80 ^ 
-4329.11 

3240.36 , 



km, = 



km. 



XI = 



-0.684779 \' 

0.476172 
, 0.55167 / 

-0.65377 \ ' 

-0.708366 

-0.266085/ 



^2 



-0.720231 \' 
-0.557622 

,-0.412703/ 
0.540953 \' 
-0.191642 

, -0.81893 / 



It is obvious that, as claimed, the separation vector between the two sites — has vanishing y-component. The 
matrix which relates the components of the vectors in the computational and earth frames is a rotation about the y 
axis: 



X^ = i?VX'^ i?V =R(0,-/3,0) =Ry(-/3) 



(6.3) 



where /3 is a time-independent (or a-independent) angle, determined by the relative orientation of the line between 
the two detector sites and the earth's axis. For the two LIGO detectors, the angle relating the Earth frame and the 
computational frame is 



/3 = -62.71383°. 

Within the computational frame, the separation vector between the two sites is 

a;?' -x^' = 2997.98km ^0^ 
and the unit-length vectors defining the arm directions are 



'0.176358\ ° 
0.476171 
,0.861486/ 



-0.69696 \ ° 
-0.557623 
0.450893 / 



X^ 



-0.536188\ ° 
-0.708366 
0.459042 / 



^2 



-0.479814\ " 

-0.191642 

-0.856185/ 



(6.4) 



(6.5) 



(6.6) 



These quantities will become useful later. 
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VII. COMPUTATION OF C 



Our goal now is to calculate Catcd as defined in (5^). To do this, we will express the spherical harmonic functions 
YifrAQ) = Yirn jO, <p) in terms of the "primed" coordinates in the computational frame. Combining the transformations 
( |6.2|) and (|6.3| ) we obtain the relationship between vectors in the cosmic and computational frames: 



=Ry(/3)R,(a)F'^ =R(a,/3,0)F^ 



(7.1) 



This transformation through Euler angles a, /?, induces a simple change in the spherical harmonics. For a given 
value of t the spherical harmonic functions in one frame are simply a sum of all the spherical harmonics with the 
same value of t in the other frame. The relation between these two sets of functions is given by the rotation matrices 



Dmm which are closely related to Clebsch- Gordon coefficients: 



(7.2) 



(See equation (4.260) of Arfkin [g[.) The inverse transformation is obtained by reversing the lower two indices on the 
rotation matrix and complex-conjugating: 



r,„(0,0)= K.fc(a,/3,o))*r^fe(^?'>'). 



The rotation matrices are conveniently expressed by equation (4.12) of Rose [Eo| 



(7.3) 



(7.4) 



Explicit formulae for the d^j, may be given either in the form of a sum, or in "summed" form. The latter expression, 
given by equation (4.14) of Rose |l^] is the most useful one for us. For m> k 



{£ - ky.{£ + m,)\ 

{e + k)i{£ - my. 



1/2 



(cosf 



2£+fe-m 



sin f 



m — k 



(m- fc)! 



2F1 {m~ e,~k~ e;m- k+1;- tan^ |). (7.5) 



Notice that because the m < i the first argument of the Gauss hypergeometric function 2-F1 is a non-positive integer 
the hypergeometric series 2F1 terminates after a finite number of terms. In fact it is possible to rewrite equation Eq. 
(7.5) in terms of Jacobi polynomials p^"-'"^*'^)^ fQj- ra>k 



(-1) 



l — rn 



/ r)\ m+k / o\ m—k 

(cosf) (-sinf) Pt:'^-''\-cos0). 



(7.6) 



1/2 / ^\ m-\-k / r)\ m—k 

{£- ky.{£ + ky._ 

In the event that m < k the d^^. may be obtained from the unitarity property, equation (4.15) of Rose 

dLkiP) = dU-f3) ^ (-l)"-'=4™(/3)- (7.7) 

Note also that the are real, so that we can drop t he c omplex conjugation that would otherwise have appeared. 

The integral over the two-sphere which appears in ( ^.8[ ) can also be expressed as an integral over all directions in 
the computational (primed) frame. In other words, J dfl — J dCl' . So our integral may be expressed as 

Cabcd{a)= ^ df„,(/?)e"""7V,'= [ sin0'd0V"-^^^^°''«'P,"(cos0') / d<j,'e''''^' QaClbn^Qd, (7.8) 
k=-e -^0 -^0 

where we have expressed the spherical harmonic functions in terms of associated Legendre functions P^, and AT = 
\xi — X2\/c denotes the light travel time between the two detector sites (10.00 msec for the two LIGO detectors). The 
normalization constants which relate the spherical harmonics and the Legendre polynomials are 



l 2i + i{i-ky 
' V An {i + ky: 



(7.9) 
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We will eventually be contracting the four indices of C with each other and with the indices of other tensors. Of course 
such contractions yield the same result in any coordinate frame, and it is easier to calculate C in the calculational 
(primed) frame. Hence we may write Cabcd{oi) = Ra"^ Rh' Rc'^ Rd^ Ca'b'c'd' where i?Q° = R2(— Q;)Ry(— /3), and 

Ca'b'c'd' ^e'""" J" di,k{f3)NH sin6iWe2"/^^™«'''p^(cos6i') / d(f>'e''"''' Cla'^ib'^ic'^id'- (7.10) 



k=~e 



The vector fi" is 



= sin 0' sin 61' . (7.11) 
V cos 61' / 



It is clear that the integral over 4>' in ( 7.1Q| ) vanishes unless k — —4, —3, • • • ,3,4. Thus even for large £ the range 



of summation over k only includes these values. There is a sense in which this reflects the fact that our signal is a 
product of the outputs of a pair of detectors, each of which has a quadrupole antenna pattern. It is also noteworthy 
that the remaining integral, over the variable 9' , can also be done explicitly for any distinct values of £ and k. 

VIII. THE REMAINING INTEGRATIONS 

We are now i n a position to evaluate the remaining integrals. We begin with the integral over a. We can re-write 
■-/im from (^.10|) as 



5 1 '•2- 



lUf) = ^7r da e-^™"e™"d?'^'(/3)4'^'(/3) x 







^a'l'c'd' dLk{P)Nt /"sin«e2"'^^^^°'"''p,^(cos0') / \4>' e'''^' (ip'tl^'^lr'^s' . 

In this integral, wc have explicitly indicated all of the dependence on a. Notice that while df^ is a function of a, 
in computational coordinates df ^ (/3) is independent of a and depends only upon (3. Likewise, the tensor Q^^'t'^c'd' 



defined by (5.9) has constant components in the computational frame. Hence the integral over a give a factor 27r: 



on 



^ti% '^rak{P)^i r sinO'de' e^^'f^^'°''>'p'',{cos0') r d(^'e^'^'hp,n,.hr'(l,,. (8.1) 

The form of this integral is interesting. The integral over 0' will vanish unless k — —4, ■ • • , 4, in which case it yields a 
product of at most four factors of sin 6*' and cos6''. Introducing a new variable u = cos 6' the two integrals appearing 
1) may be expressed as linear combinations of the integrals 



' du e2-/^^«P,^(u)u^(l - u^p/^ (8.2) 
-1 

where is a non-negative integer bounded by -1- |fc| < 4. Such integrals can be expressed in closed-form as we show 
in the Appendix. 

From this point on, wc consider only the case of the two LIGO detectors. In this case the simplest way to proceed 
is to compute 

sk{e') ^ df\p)d^^'i\p)ei,i,l'i r dcb'e'''^'Cip,ng,Cir'(is'. (8.3) 

Jo 

To evaluate this, we use definition (|4.5D of the df'', the contraction operator (^.9D and the arm directions (|6^). 



Substituting in the vector D,"" given by (7.11) gives elementary integrals over 0'. The results are easily written in 
terms of the variable u = cos 6': 
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si{u 

S2iu 
S3{U 
S4{u 
S-k{u 
Sk[u 



= -3.01308 + 1.75421 + 0.945109 

= (0.5451 + 0.543353 i - (1.6954 + 1.73284 i) u^) u{l~ u^f/^ 
= (-0.0428705 + 1.41125 i + (0.0119604 - 0.670812 i) u^) (1 - u^) 
= (-0.245744 + 0.227815i) u (1 - u^f/"^ 
= (0.0492709 + 0.00345516i) (1 - u^ f 

= for \k\ > 4 



Wc then have 



(8.4) 



We now evaluate these functions for the first few muhipoles. For this purpose we introduce a dimensionless frequency 
variable x = 2nfAT. Because = (— 1)^~''™7^ we give these functions only for m = 0, • • • , ^. They may be 
conveniently written in terms of spherical Bessel functions j„. For £ = one has 

70 o(a;) = -0.0352174jo(a;) - 0.818115^^^^ + 0.848647^^^. 



A graph of this function is shown in Figure 

Y 



Freq f, Hz 



300 




FIG. 1. The (real) function 70, o(/) is shown for the LIGO pair of detectors. 



For 



1 one has 



71 o(a;) = -0.0279637ui(a;) - 0.252119^^^^ - 1.66955?:^^ 

X x^ 

7i,i(a;) = 0.0383329iji(a;) - (0.327033- 1.03547i):^^ + (1.90568 - 1.77847z)^^. 



A graph of these functions is shown in Figure | 
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FIG. 2. The functions 71, m(/) are shown for the LIGO pair of detectors. The real parts are shown as the solid curves, and 
the imaginary parts as the dotted curves. 
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For 



2 one has 



72,o(a;) = 0.0145494jo(a;) + 1.00009 



9.39901^^^7^ + 28.0344:^Ar^ 



72,1 (x) = 0.0392947jo(a;) + (0.385015 + 0.335238i) 



-(6.01443- 16. 1191i 



72,2(2;) = -0.0380887jo(x) - (1.05867 + 0.649899i) 
-(18.5056 + 20.4538i):^^. 



ii(a^) 



(2.38288 + 2.485340^^^ 



(9.05483 + 6.80403i):^^ 



A graph of these functions is shown in Figure | 
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FIG. 3. The functions 72, m(/) are shown for the LIGO pair of detectors. The real parts are shown as the solid curves, and 
the imaginary parts as the dotted curves. 

For £ = 3 one has 

73 o(a;) = 0.0416301iji(a;) + 0.805209^:^^ - 11.2779^:^^ + 18.399i^i^ 

X x'^ x-^ 

73 i(x) = 0.00182185ui(x) - (0.0155429 + 1.03107i):^^ - (1.38339 - 11.1549^):^^ 

X x'^ 

-(7.72558 + 45.6531^):^^ 

x-^ 

73 2[x) = -0.0461979iji(a;) + (0.788263 - 0.534988^):^^ - (7.13839 - 8.50643^):^^ 



+ (42.9458- 1.80521i):^^ 

X'^ 

73 3(2;) = 0.0365627iji(a;) - (0.935792 - 1.08142^; 
-(33.6244- 32.4339i):^^. 



+ (11.3339 - 13.8287i):^^ 



A graph of these functions is shown in Figure 
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FIG. 4. The functions 73,m(/) are shown for the LIGO pair of detectors. The real parts are shown as the solid curves, and 
the imaginary parts as the dotted curves. 

For f = 4 one has 

74 n(a;) = 0.0232299jo(a;) - 0.899008^^1^ + 9.78424^^^ - 88.696:^^ + 438. 555^^1^ 



74^i(a;) = -0.0367891jo(a;) - (1.10932 + 0.313862i 



(27.3542 + 5. 98172i):^^ 



-(204.205 + 48.3679i):^^ + (496.286 + 149.558?):^^ 

x"^ x^ 

74 2(a;) = -0.0155407jo(a;) + (1.05639 - 0.265167^):^^^ - (14.276 - 0.635514^):^^ 

X x^ 

+ (127.588 - 28.9739i)^^^ - (551.647 - 310.656i):^^ 

x^ 

74 3(2;) = 0.0502848jo(a;) + (0.521081 + 1.287j)^^^ - (20.3464 + 20.2033i)^^^ 

X x^ 

+ (135.553 + 177.464^):^^ - (156.429 + 719.09H):^^ 

X'^ x^ 

74 4(2;) = ~0.0344655jo(a;) - (0.977416 + 1.17615i):^i^ + (24.322 + 22.0573^):^^ 

-(183.738 + 163.238z):^^ + (458.028 + 438.49i)^^^ 

x° x^ 

Note that the jim with £ odd vanish as / ^ 0, in contrast with the functions with £ even, which approach constant 
values at zero frequency. 



IX. OBSERVABILITY AND SIGNAL-TO-NOISE RATIOS 



Up to this point, we have shown how anisotropics in the gravitational- wave stochastic background give rise to 
periodic variations in the "signal" obtained by correlating a pair of detectors. These periodic variations are described 
by the Fourier series (5.1), with coefficients Sm- In this section, we address the question: how precisely can the 
values of these coefficients be determined, in the presence of noise in the two detectors? We answer this question by 
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calculating the signal-to- noise (S/N) ratios that would arise in measurements of the Sm, which also permits us to 
determine the best choice of the optimal filter functions Q{f). 

In carrying out this analysis, we follow the technique used in Section 3.2 of reference ||^. The following presentation 
will be somewhat crypic as we will assume that the reader is familiar with that material. 

The noise in the detectors is characterized by a cross-correlation function 

{nKf)nAf')) = ls.Af-nmf\)- (9.1) 

Here i — 1,2 labels the two detector sites and Pi{f) is the (one-sided, real) noise power spectrum of the I'th detector. 
For the initial and advanced LIGO detectors, these power spectra are shown in Fig. ^. For our calculations of S/N 
we will need to know the noise properties of the detectors averaged over our "windowing time" r. To obtain these, 
we first characterize the detector noise in the time domain, by Fourier transforming ( |9.1| ) . This gives 

/OO pOQ 
dfe-'^^f' / d/'e2"/'*'«(/)n,(/')) 
-OO J —OO 

1 z*^ 

-2^vJJfe-'-'^^'-'^m\f\)- (9.2) 
In words, this says that the Fourier transform of the noise auto-correlation function is the noise power spectrum. 
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FIG. 5. The predicted noise power spectra of the initial and advanced LIGO detectors. The horizontal axis is logio of 
frequency /, in Hz. The vertical axis shows logio(P(f)/sec)^'^^, or strain per root Hz. These noise power spectra are the 
published design goals. The bumps appearing in the low-frequency part of the advanced LIGO noise curve are obtained by 
folding measured seismic noise data with the predicted transfer function of the seismic isolation (stack) system. 



We can now find the "windowed" version of these formulae, using the definition (4.6) of the windowed transform. 
Taking the windowed transforms of (9.2) yields 



(n*(/,i)n,(/',i')) 



t+T/2 



t'+T/2 



dt" e^'^'f'" / dt'" {n,{t")n,it"')) 

t-T/2 Jt'-T/2 

i5„e2"(/*-/'*') 
2 ^ 



dpe-^-^P^'-''^5Af - p)Sr{f - P)m\p\)- 



(9.3) 



Note that in the long averaging time limit r ^ oo this reproduces equation (^T|). 

In virtually any realistic scenario, the intrinsic instrumental detector noise is expected to be much larger than the 
strain arising from the stochastic background of gravity waves. For this reason, if we define the "noise" in a given 
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measurement of Sm by Nm = Sm — (Sm) then to g ood approximation a formula for Nm may be obtained by replacing 
the total detector output Si{t) which appears in (42) by ni{t). Using the definition of the "signal" (49) and the 
definition of the m'th harmonic (j^^ one obtains 

1 f°° 

Nm^Tf dte~'' dfhl{f,t)n2{f,t)Q{f). (9.4) 

^ Jo J-oc 

(Note that in this formula, we assume that the total observation time T is large compared with the earth's rotation 
time Te and that during the observation time the earth has made an integral number of rotations, so that T/Tc is a large 
integer. We will typically take T to be 1 sidereal year.) The approximation that we make here is obviously consistent 
with (Nm) — since the noise in the two detectors is assumed to be uncorrelated, so that {nl{f,t)h2{f',t')) = 0. 

The "noise" arising in the measurement of Sm may be characterized by (|iVmP) = (|<5'mp) — K^*™)!^. We now 
calculate {N*^Nm')- Substituting equation (9.3) into equation (9.4) and "squaring" yields 

{N*^Nm')^Tf^ dt rft'e-eM-™'0 / dfQ*{f) dfQ{n{nlU,t)Mr,t')){nlU,t)fi2{rX)) 

^ Jo Jo J-oo J-OG 

1 pT pT pOO pOQ pOO pOQ 

= — 2/ / di'/ df df / dp dp'Q*(/)Q(/')e*'^'=("*-™'*')+2"(p-p')(t-t') 

Jo Ja J-oo J-oo J-oo J-oo 

X Srif - p)Srif' - p)Srif - p')5r{f' - p')Pl{\p\)P2{\p' \) ■ (9.5) 

We can simplify this expression to the point where it is useful, however this is somewhat tricky - there is only one 
order in which the integrals above can be simply evaluated to yield useful approximations. We first do the integrals 
with respect to t and t' exactly. This gives 

1 j*oo poo />oc poo 

{NmNm') = -^ / df / df / dp / dp'g*(/)0(m(b|)P2(b'|)e--"^("-™')/^ 

J —OO J —oo J —oo J —oo 

X 6t{p -p' + mu;j27r)dT{p - p' + m'uj2^)5r{f - p)5r{f' - p)5r{f ~ p')5r{f' ~ p')- (9.6) 

Now we note that the effective support of Srif) extends over a very narrow range of frequencies (typically, |/| < 10~^ 
Hz) compared with the effective support of 5r{f) (whose support is typically |/| < 10~^ Hz). In addition, none of the 
remaining integrand varies over such a narrow frequency range. So we are justified in replacing St^P ^ p' + muj^ /2tt) 
by the ordinary Dirac delta function 5{p — p' + mti>o/2vr). This gives 

-1 pOO POO POO 

(iV;;iV™,) = — 2 / df df Q*{f)Q{f') dp'<5,(/-p')'5r(/'-p')e-^"°^^™-'"')/' 



LUc{m ~ ■m')\ „ /. , mciJc\ ^ ( j-i i rnujc 



2n ' V 2tt J V ^ 27r 



^Sr( ^'''':''^' ]Sr(f-p' + ^]Sr(f-p' + ^]P^(p'-^ ]PM). (9.7) 



2tt 



If we now additionally assume that the observation time T is much greater than the period of a single rotation — ^ 
then to good approximation St ^ ""("^"^ ) ^ ~ Tdmm'- Thus 

-t pOO POO POO 

(iV;„iV„,) = — <5w / df / df Q*{f)Q{f') / dp Srif - P)W - P) 

J —oo J —OQ J — oo 

^Sr{f-p+^)sr{r-p+ ^) p.{\p~^\) mp\)- (9.8) 

Next we note that the width of the Sr (/) in frequency space is quite large compared to mujc provided that we restrict 
m to be fairly small: |m| < To/r. We also assume that the noise power spectrum Pi{f) does not vary significantly 
over frequency scales of lOc ~ 10^^ Hz. In this case we can neglect the shifting of the arguments by mujc above, 
obtaining 

{n:^N„,) = --<5w / df / df Q*{f)Q{f')SUf - f')Pi{\f\)P2{f\) 

^ ^ J —oo J —oo 

t2 /"^ 

^Smm'J^ d/|Q(/)pPi(|/|)P2(|/|). (9.9) 

J —oo 

Setting m — m' we finally obtain an expression for the expected "squared noise" in a measurement of Sm- 
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{\Nm\') 



T 

4T 



df\Q{f)fp,{\f\)P2m- 



(9.10) 



We can make use of this expression to find the optimal filter Q{f). 

The (squared) Signal-to-Noise ratio S/N for the m'th harmonic is no w d efined via the ratio of expect ed si gnal 
(magnitude squared) divided by expected (squared) noise. Making use of (5.3) for the former quantity, and ( 9.1C| ) for 
the latter yields 



4T| 



T C^dfQ{f)H{f) J2 Ptmlimif) 

e=\m\ 



•r^dfmwp.iifmm 



(9.11) 



Notice that the averaging time r (which was earlier chosen in a rather arbitrary manner) drops out of this expression. 
Provided that the assumptions about r (|4.l|) used in deriving this equation are satisfied, the actual value is irrelevant. 

In order to find the optimal filter function, it is useful to introduce an inner product. For any complex functions of 
frequency A{f) and B{f), this defines a complex number which is denoted by {A, B). The definition is 



/oo 
dfA*{f)B{f)p,{\f\)P2m- 
-oo 



(9.12) 



This inner product is positive definite because {A, A) is real and non-negative, vanishing only if A is zero. In terms 
of this inner product, the Signal-to-Noise ratio may be written as 



S 



Svr 



'5' Pi{\f\)P2{\f\) X! P*emlemif) 

l=\m\ 



iQ,Q) 



(9.13) 



The optimal choice of filter function Q{f) for determining the m'th harmonic is the one which maximizes this ratio. 
The largest value is obtained by choosing 



H{f) 



Pli\f\)P2{\f\) 



(9.14) 



f=|m| 



Using the definition ( |3.11 ) of H{f) in terms of the spectral function r2gw(/), and substituting the optimal filter Q.^, 
into the expression for S/N yields 



9^ 
507r4 



T 



df 



PPlif)P2if) 



^ Plnijeniif) 
?=|m| 



(9.15) 



For any given source of stochastic gravitational waves, one can use this formula for {S/N)m to determine the obser- 
vation time T required to observe the m'th harmonic of the signal as the earth rotates relative to the cosmological 
frame. 

In precise analogy with the analysis given in the m'th harmonic is observable with 90% confidence if {S/N)m 
exceeds 1.65. 



X. EXAMPLE: DIPOLE INDUCED BY PROPER MOTION 



It is well known that the electromagnetic background radiation, generally referred to as the Cosmic Microwave 
Background Radiation (CMBR), is highly isotropic. The largest deviation from isotropy results from the motion of 
our local system (the solar system barycenter) with respect to the cosmological rest frame. Analysis of data from the 
Cosmic Background Explorer (COBE) satellite shows that our local system is moving with a velocity /3propcr = v/c = 
0.001236 in the direction (I = 264°, 6 = 48°) in galactic coordinates, or equivalently (a = 168°, (5 = —7°) in celestial 
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coordinates. To lowest order in our proper velocity, this gives rise to an anisotropy in the CMBR described by the 
temperature distribution 



T(7) =To(l+/3properCOS7), 



(10.1) 



where 7 is the angle between a point on the sky and the velocity vector of our local system, and Tq = 2.73K is the 
mean temperature of the CMBR. 

In this section, we address the question: "Would a corresponding dipole moment in the stochastic gravity wave 
background, arising from the proper motion of the solar system barycenter, be observable with either the initial or 
advanced LIGO detectors?" In the previous section, we calculated the Signal-to-Noise ratio for observations of the 
fn til harmonic Sm (with respect to the earth's rotation) of the signal (obtained by correlating two gravitational-wave 
detectors). This is determined entirely by the quantities 

Hq = Hubble expansion rate, scc~^. 
T = Total observation time, sec. 
Itmif ) = Overlap reduction functions, dimensionless. 

f2gw = Fractional energy density in gravitational waves, dimensionless. 
Pim = Normalized multipoles of gw stochastic background, dimensionless. 
Pi = Noise power spectral density of detector i, sec. 

As an example, we use this formula to answer the question posed above. 
Let us make the following reasonable assumptions: 

(1) rig„(/) is constant in the LIGO band, and 

(2) the stochastic gravitational-wave backgroimd is isotropic in the same rest frame as the CMBR. 
The anisotropies in the stochastic background resulting from our proper motion are then described by 



PemXtm (^, (A) = 1 + /^proper (cos 6 COS 97° + sin 6 sin 97° cos (j)) 

im 



(10.2) 



where the angles above are standard spherical coordinates in the (barred) earth frame. One thus obtains the following 
multipole moments 

Poo 
Pi -1 



/3proper\/^sin97 



Plfi = /JproperV — COS 97° 



Pl,l 
pern 



= -/3p 




sin 97° 



for £ > 1. 

To detect this signal, the optimal filter functions are 

^"^•^^=^™/3Pl(|/|)P2(|/|) 



(10.3) 



(10.4) 



e=\m\ 



where Cm is an (irrelevant) normalization constant. Making use of the optimal filter functions Qq and Qi for the 
monopole and dipole terms, we can make predictions about how large figw needs to be in order that Sq and Si are 
observable with a given level of confidence in a given observation time. 

Wc can express the sensitivity of a search for the m'th harmonic in terms of the minimum value of fig^ necessary 
to observe it with 90% confidence. For 90% confidence we need a Signal-to-Noise of 1.65. The minimum value of figw 
is then given for the m'th harmonic by 



df 



^ PemHm.{f) 

e.=\m\ 



(10.5) 
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These values are shown in Table |y for the initial and advanced LIGO detectors. We note that there are potential 
sources of sufficient intensity that a dipole might be observable with the advanced LIGO detector. These include 
stochastic backgrounds due to cosmic strings or to a population of unresolved cosmological-distance supernovae. 



TABLE I. Sensitivity of Initial and Advanced LIGO detectors to a dipole term in the gravitational stochastic background, 
arising from motion of our local system. This table shows the intensity of stochastic background required to detect either the 
monopole (So) or dipole [Si) term in the signal, with 90% confidence, in one year of observation. 





Initial LIGO 


Advanced LIGO 


Monopole (m=0) 


ngo% = 1.6 X 10-^/1^00 


^90% = 1.6 X 10-"/i-oo 


Dipole (m=l) 


ngo% = 2.5 X 10-3/1^00 


^90% = 5.3 X 10-»/iroo 



XI. FILTER NORMALIZATIONS AND SIGNAL STRENGTHS 



Specific models for an anisotropic background predict signal harmonics Sm of definite amplitudes and phases. The 
phase of a given Sm depends upon the chosen origin of time; shifting the origin of time changes the phase of Sm but 
not its amplitude. In order to falsify or confirm a particular anisotropic model, one needs to predict the set of complex 
numbers Sm- However these numbers depend upon the normalization Cm of the optimal filter function Qm{f) so in 
order to make predictions, we need to adopt a convention for these. 

We have already chosen optimal filters for which the expected value of Sm is real. At least upon first inspection, 
there do not appear to be any very convenient choices for the overall scale of normalization. Here, for illustrative 
purposes we adopt (what appears to us as) the least disagreeable of these. We choose the normalization Cm so that 
the expected amplitude of the noise^ is unity: {\Nm\'^) = 1- We also assume that C,„ is positive and real. Note that 



with this choice, the "signal" is dimensionless. From equation (9.1C) this gives 



Cm. 



2Vt 



df 



pm)P2i\f\) 



t=\m\ 



(11.1) 



Having chosen a normalization, we can list the expected signals for a given detector. If the actual signal value exceeds 
1.65 then it has been detected wit h 90% confidence. For example, in the case of the dipole source just discussed, the 
expected signals are given by (5.3) as 



{S„ 



df 



-00^(1/1)^2(1/1) 



Evaluating this for our previous example gives: 



Ptynllraif) 



(11.2) 



1.65, 



1.6E; 



1 year 1.6 x 10^^ 



T 



Jigw'tioo 



1 year 1.6 x lO"" 



for initial LIGO 
for advanced LIGO 



1.65, 



1.65, 



1 year 2.5 x lO^^ 

T flg^h^QQ 



1 year 5.3 x 10" 
(Sm) for TO > 2 



for initial LIGO 
for advanced LIGO 



(11.3) 
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XII. GALACTIC SOURCES 



In this section we consider the possibihty of detecting anisotropics in the stochastic gravitational wave background 
assuming that this background either originates in, or is scattered in the same way as, the luminous matter in 
our galaxy, and for sources distributed in the same way as the galactic halo. It appears very unlikely that in the 
LIGO/VIRGO/GEO frequency band there are any sources of a stochastic gravitational- wave background distributed 
in this way. However the absence of such harmonics would be one way to demonstrate that a stochastic background 
had extra-galactic origin. 

For the distribution of luminous matter in the galaxy we consider a set of three models constructed by Kent, Dame 
and Fazio [|l^ to model 2-4 /im data from the Infrared Telescope taken as part of the Spacelab 2 mission. These models 
all assume cylindrical symmetry and take the total luminosity to consist of two components, one, i/£)(r, z), modeling 
the disk and the other, vsir, z), modeling the central bulge, where r and z denote cylindrical polar coordinates based 
at the center of the Galaxy. 

The first model, which we refer to as KDFl, takes 

sech^(z/(2fe,)) 

lyoir, z) ^ ^ioe ^— (12.1) 



o{[r' + {z/{l-eB))Y'/hB) 



iyBir,z)^fiB ^ r (12-2) 

nhB 

where fio = 1072 Lgpc-^ K = 2775 pc, = 121pc, ^J.B = 6208 LqPC"^^ Hb = 634 pc and eb = 0.26. Here Kq 
denotes the modified Bessel function or order 0. 
The second model (KDF2) takes 

exp(-|z/(2/i,)|) 

,^D{r, z) - Aice-'^/"- ^ (12.3) 



K,{\r^ + {z/{l-eB)tf'/hB) 



VB{r,z)=^iB ^ r (12-4) 

where /in = 1208 Lopc-^, hr = 2694 pc, = 204 pc, = 771OL0pc-2, Kb = 500 pc and es = 0.19. 
The third model (KDF3) takes 



cxp \—\z 



/(2/i.(r))|) 



VD{r.z) = ^iDe-^"'^■ ^ ^, . . ^ (12.5) 

2/i2(r) 



K 



.{[r' + {zl{l-eB)f]"'/hB) 



i^B{r,z)=fiB r (12.6) 

nhB 

where fijj = 978Lqpc"^, hr — 3001 pc, fiB = 7395 Lqpc"^, hB — 667 pc, eb — 0.39 and /iz(r) is now a function given 
by 

/imin + (/l0 - /lmin)7 T r > 



' min 



with hjnin = 165 pc, Tijiin = 5300 pc, hQ = 247 pc and rQ = 8000 pc, the last two quantities being the height of the 
galactic disk at the position of the Sun and the distance of the Sun from the Galactic center, respectively. 
Finally, as a model for the galactic halo we take the model of Young 

exp [-7.669(i?/i?e)^/^] 

I'halo = Mhalo / „ , „ sy/o R > 0.2i?e (12-8) 

(U/Ue)''° 

where R = {r^ + z^Y^^ denotes the distance from the Galactic center and Re = 2700 pc. 



19 



The final piece of information we need before we can calculate the multipole moments pim for such a distribution 
is the direction to the Galactic center and the orientation of the Galaxy. It is standard to express this in terms of the 
equatorial coordinate system in which the z-axis is taken along the Celestial North pole, the x-axis is taken in 
the direction of the (1950) vernal equinox. 90° — 6 and a are taken as the spherical polar coordinates corresponding 
to these axes. In terms of these, the direction to the Galactic center is given by 

a = 265.6° 5 = -28.9°, (12.9) 

and the direction of the Galactic North pole by 

a = 192.25° (5 = 27.4°. (12.10) 

We take the equatorial coordinate system to define our Cosmic frame. We now determine -P(f2) by choosing a 
direction (l and integrating the intensity along that direction. The number of sources in an element of solid angle 
dil at distance is D from the Earth is proportional to v D^dDdQ.. The intensity of the source drops like l/Z?^ and 
hence t he in tegrated intensity is just proportional to / vdD. In this way we arrive at pim for each model and hence 



by Eq. (10.5) the minimum value of fJgw necessary for detection of each multipole. To ensure convergence of the sum 



in this case it was necessary to include contributions to the sum for I up to of order 100 (dependent on m). Tables |^ 
and III show the intensity of the stochastic background distributed in this way required to detect multipole moments 
from m = to m = 24 for the four models with 90% confidence in one year of observation for Initial and Advanced 
LIGO respectively. 



TABLE II. Sensitivity of the Initial LIGO detector to the first 25 multipoles in the gravitational stochastic background, 
assumed to follow the luminosity of the Galaxy or Galactic halo. This table shows the intensity of stochastic background 
required to detect the multipole Sm with 90% confidence in one year of observation. 

m h\ooQ,gQ% KDFl /i?oo^^90% KDF2 /i?oof^90% KDF3 /i?oof^90% Halo 






1.4 


X 


10^5 


1.4 


X 


10-5 


1.4 


X 


10-5 


3.9 


X 10-^ 


1 


1.4 


X 


10"5 


1.3 


X 


10-5 


1.3 


X 


10-5 


3.4 


X 10-^ 


2 


1.9 


X 


10"5 


1.9 


X 


10-5 


1.8 


X 


10-5 


3.6 


X 10-^ 


3 


3.0 


X 


10"5 


2.9 


X 


10-5 


2.9 


X 


10-5 


4.2 


X 10-^ 


4 


4.6 


X 


10~5 


4.4 


X 


10-5 


4.6 


X 


10-5 


5.0 


X 10-^ 


5 


7.3 


X 


10~5 


6.9 


X 


10-5 


7.1 


X 


10-5 


6.1 


X 10-^ 


6 


7.8 


X 


10~5 


7.4 


X 


10-5 


7.5 


X 


10-5 


8.1 


X 10-^ 


7 


8.6 


X 


10"5 


8.0 


X 


10-5 


8.3 


X 


10-5 


1.2 


X 10-^ 


8 


1.2 


X 


10-'' 


1.1 


X 


10-* 


1.1 


X 


10-* 


1.7 


X 10-^ 


9 


2.0 


X 


10-'' 


2.0 


X 


10-* 


1.9 


X 


10-* 


2.7 


X 10-^ 


10 


3.6 


X 


lo--* 


4.0 


X 


10-* 


3.5 


X 


10-* 


4.3 


X 10-^ 


11 


5.0 


X 


10-* 


5.2 


X 


10-* 


4.9 


X 


10-* 


7.0 


X 10-^ 


12 


6.3 


X 


10-* 


6.2 


X 


10-* 


6.2 


X 


10-* 


1.2 


X 10-^ 


13 


9.2 


X 


10-* 


8.9 


X 


10-* 


9.0 


X 


10-* 


2.0 


X 10-^ 


14 


1.5 


X 


10-^ 


1.5 


X 


10-^ 


1.5 


X 


10-^ 


3.1 


X 10-^ 


15 


2.4 


X 


10-^ 


2.5 


X 


10-^ 


2.4 


X 


10-^ 


5.0 


X 10-^ 


16 


3.6 


X 


10-^ 


3.7 


X 


10-^ 


3.6 


X 


IQ-^ 


7.8 


X 10-^ 


17 


5.1 


X 


IQ-^ 


5.2 


X 


10-^ 


5.1 


X 


IQ-^ 


1.2 X 10" 


18 


7.3 


X 


IQ-^ 


7.2 


X 


10-^ 


7.1 


X 


10-^ 


1.7 X 10" 


19 


1.0 


X 


10-^ 


1.0 


X 


10-^ 


1.0 


X 


10-^ 


2.4 X 10" 


20 


1.4 


X 


10-^ 


1.4 


X 


10-^ 


1.4 


X 


10-^ 


3.6 X 10" 


21 


2.1 


X 


10-^ 


2.1 


X 


10-^ 


2.0 


X 


10-^ 


5.0 X 10" 


22 


2.9 


X 


10-^ 


2.9 


X 


IQ-^ 


2.8 


X 


IQ-^ 


6.5 X 10" 


23 


4.0 


X 


IQ-^ 


3.9 


X 


10-^ 


3.8 


X 


IQ-^ 


8.9 X 10" 


24 


5.1 


X 


10-^ 


5.0 


X 


10-^ 


4.9 


X 


IQ-^ 


1.2 


X 10^ 
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TABLE III. Sensitivity of the Advanced LIGO detector to the first 25 multipoles in the gravitational stochastic background, 

assumed to follow the luminosity of the Galaxy or Galactic halo. This table shows the intensity of stochastic background 



required to detect the mu 


Itipole Sm 


with 90% confidence 


in one year 


of observation. 










m 


hioo^go% KDFl 


/if 00^1 


,0% KDF2 




)o% KDF3 


^100^90% Halo 





1.8 X 10-1° 


1.7 X 10-1° 


1.7 X 10-1° 


6.7 X 


10-® 


1 


4.8 X 10-1° 


4.7 X 10-1° 


4.7 X 10-1° 


7.8 X 


10-® 


2 


1.1 X 


10-3 


1.0 


X 


lo-'' 


1.0 


X 


10-8 


1.1 X 


10-^ 


3 


5.3 X 


10-^ 


5.0 


X 


io-« 


5.4 


X 


10-8 


2.6 X 


10-^ 


4 


6.4 X 




6.0 


X 


io-'> 


6.3 


X 


10-8 


6.4 X 


10-^ 


5 


1.1 X 


10-® 


1.0 


X 


10-® 


1.0 


X 


10-® 


1.6 X 


10-® 


6 


3.1 X 


10-* 


3.0 


X 


10-® 


3.0 


X 


10-® 


4.4 X 


10-® 


7 


1.2 X 


10-^ 


1.3 


X 


10-^ 


1.2 


X 


10-^ 


1.3 X 


IQ-® 


8 


2.4 X 


10-^ 


2.4 


X 


10-^ 


2.4 


X 


10-^ 


2.7 X 


10-® 


9 


4.5 X 


10-^ 


4.2 


X 


10-^ 


4.5 


X 


10-^ 


6.3 X 


10-® 


10 


1.1 X 


10-® 


1.0 


X 


10-® 


1.1 


X 


10-® 


1.3 X 


10-* 


11 


3.0 X 


10-® 


3.3 


X 


10-® 


3.1 


X 


10-® 


2.6 X 


10-* 


12 


5.5 X 


10-® 


5.5 


X 


10-® 


5.7 


X 


10-® 


4.5 X 


IQ-* 


13 


8.3 X 


10-® 


7.7 


X 


10-® 


8.7 


X 


10-® 


7.7 X 


10-* 


14 


1.5 X 


10-^ 


1.5 


X 


10-^ 


1.6 


X 


10-^ 


1.3 X 


10-^ 


15 


3.0 X 


10-^ 


3.2 


X 


lo-'' 


3.1 


X 


10-^ 


2.0 X 


10-^ 


16 


4.4 X 


10-^ 


4.4 


X 


lO"'^ 


4.7 


X 


10-® 


3.0 X 


IQ-^ 


17 


6.0 X 


lO"""' 


5.6 


X 


10-'' 


6.4 


X 


10-® 


4.4 X 


IQ-^ 


18 


9.4 X 


10-^ 


9.0 


X 


10-^ 


1.0 


X 


10-* 


6.4 X 


10-^ 


19 


1.6 X 


10-" 


1.6 


X 


10-* 


1.7 


X 


10-* 


9.0 X 


10-^ 


20 


2.1 X 


10-* 


2.1 


X 


10-* 


2.3 


X 


10-* 


1.2 X 


10-2 


21 


2.7 X 


lo-'' 


2.5 


X 


10-* 


3.0 


X 


10-* 


1.7 X 


10-2 


22 


3.9 X 


10-* 


3.7 


X 


10-* 


4.2 


X 


10-* 


2.3 X 


10-2 


23 


5.9 X 


10-* 


5.8 


X 


10-* 


6.4 


X 


10-* 


3.0 X 


10-2 


24 


7.8 X 


10-* 


7.5 


X 


10-* 


8.4 


X 


10-* 


3.9 X 


10-2 



XIII. CONCLUSION 

In this paper, we have shown how the signals from a pair of gravitational wave detectors may be analyzed to search 
for anisotropies in the stochastic gravitational wave background. We have shown how the correlation between two 
detectors may be determined with an averaging time short compared to a day but long compared to the light travel 
time between the detectors, and how this correlation may be decomposed into harmonics of the Earth's rotation. 
We have calculated the signal-to-noise ratios associated with such measurements, and shown that certain types of 
anisotropy might reasonably be detected with instruments that will be available in the not-too-distant future. For 
example the anisotropy due to our motion with respect to the cosmological rest frame might well be observable with 
a second-generation LIGO detector. 

At this point, it is difficult to proceed further without a more detailed knowledge of the instrumental data that will 
be forthcoming. The results given in this paper make it straightforward to predict the expected harmonic amplitudes 
Sm of the detector correlation for a given anisotropic distribution p£m of gravitational wave background. It is more 
difficult to go the other way. This would involve using a given set of observed harmonic amplitudes Sm to obtain 
the values of (or constraints on) the pim- There are a variety of fitting techniques that could be used - making the 
appropriate choice will probably require real data. 
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APPENDIX: ANALYTIC EVALUATION OF INTEGRALS 

In this appendix we derive a closed form expression for the integral 

Ie{x)=J dMe*""Pf(w)w^(l-M2)|fc|/2 (A^) 

where iV is a non-negative integer. 
First we note that 

PJ^iu) = |_||:p,^(^,) (A2) 

so that it is only necessary to deal with the case fc > 0. For fc > we have 

pUu) = (-l)'^(l - u^f'^^P,{u) (A3) 

I'lix) = due'^-u''{l-u'f^P,{u). (A4) 

The presence of the factor {1 — v?')^ now ensures that when we integrate by parts k times no boundary terms appear 
so that 

I^x) = du Pi{u)^ [e™^M^(l - u^f] . (A5) 
The derivative can be expanded by the Liebniz rule to give 

k 

Y,{ixrpN+k+r{u)e'''^ (A6) 



givmg 



r=0 



where pN+k+r{u) is a polynomial of degree N + k + r. Our problem is thus reduced to that of finding 

J^\x)^ j due'^'^'u^^Piiu). (A7) 

For AI — this is an elementary integral given by 

J,"{x) = 2i'jfXx) (A8) 

where je{x) denotes the spherical Bessel function of order £ jl^ . Higher M values may then be obtained by differen- 
tiation 

These derivatives may in turn be expressed back in terms of (undifferentiated) spherical Bessel functions using the 
relations ||l^ 

^M^) = W^-ii^) l)je+i{x)] (AlO) 

i 

= -Je{x)-Ji+i{x). (All) 
x 

For M = 1 to 8 we have 
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{2£+iy''~'^ ' (2£+l)' 

{2e + 21-1) + + 2) . 



(2^-l)(2^+l)-''"^' ' (2^- l)(2^ + 3)''*' ' (2£+l)(2£ + 3)" 



(2£ - 3)(2£ - 1)(2£ + 1)-"^-^^ ' (21 - 3)(2£ + + 3) 

3(£+l)(£V2£-l) . (^+l)(£ + 2)(£ + 3) . 

(2^ - 1) (2f + 1) (2^ + 5) •^^+^ (2^+1) (2^ + 3) (2f + 5) 

(.{l-X){l-2){l--i) , , 2^(^-l)(2£2_2^-7)) 



(2£ + - - 3)(2£ - 5)-'* ' ' {21 + 3)(2^ + \)[2i - 1)(2£ - 5)" 
3(2^4 _^ 4^3 - 6^2 _ _^ 3) 2(£ + l)(£ + 2)(2£2 + 6^-3) . 



(2£+5)(2£ + 3)(2£- l)(2£-3)-' ' ' (2£ + 7)(2£ + 3)(2£ + 1)(2£ - 1) 
(^+l)(^ + 2)(^ + 3)(f + 4) 



(2f + 7)(2£ + 5)(2£ + 3) (2^ + 1) 



l)(£-2)(£-3)(^-4) . 5^(£- l)(£-2)(£2 -2^-5) . 



(2f-7)(2f-5)(2^-3)(2^-l)(2£ ^ (2^ - 7)(2£ - 3)(2£ - 1)(2£ + 1)(2£ + 3)" 
5^(2^4-16^2 + 23) , , 5(^+l)(2^4 + 8^3_4^2_24^ + 9) 



(2£-5)(2£-3)(2£+l)(2^ + 3)(2^ + 5)-' ' ' (2£ - 3)(2^ - 1)(2£ + 1)(2£ + 5)(2^ + 7) 
5(^+l)(^ + 2)(^ + 3)(^2^4^~2) . , , (^+l)(£ + 2)(£ + 3)(£ + 4)(£ + 5) 



(2f - 9)(2f - 7)(2f - 5)(2£ - 3)(2£ - l)(2f + 1) 

3£(£ - 1)(^ - 2){l - Z)(2e -U- 13) 
{21 - 9)(2£ - 5)(2£ - 3)(2^ - 1)(2^ + 1)(2£ + 3) 

\U{1 - \){t - 2^ - we + 12£ + 29) 
(2^- 7)(2^- 5)(2^- 1)(2£+ l)(2£ + 3)(2^ + 5) 
5(4£6 + 12^5 _ 3g^4 _ gg^3 ^ gg^2 ^ ;^3g^ _ 

{21 - 5) {21 - 3) (2£ - 1) (2^ + 3) {21. + 5) (2^ + 7) 

\h{l + 1)(^ + 2){t + + - 2U + 9) 
(2£ - 3) (2£ - 1) (2£ + 1) (2£ + 3) (2i? + 7) (2£ + 9) 

3(^ + 1)(^ + 2){l + 3)(£ + 4)(2£2 + 10^ - 5) 
{21 - \){2l + 1)(2£ + 3)(2£ + h){2i + 7)(2£ + 11) 

{i + \){l + 2)(£ + 3)(£ + 4)(£ + 5)(£ + 6) 
{21 + l)(2i' + 3)(2^ + 5)(2^ + 7)(2^ + 9)(2^ + 11) 

/((-i)((-2)i(-3)(r-4)(r-5)(/-(.i) 

(2^- ll)(2£-9)(2^-7)(2£-5)(2£-3)(2f- 1)(2^+1) 
lt{l-\){l-2){l-Z){t-\){e 



3e-6{x) 



(2£- 1)(2£+ l)(2<' + 3)(2^ + 5)(2^ + 9)-' "^^ ^ (2£ + 1)(2£ + 3)(2^ + 5)(2£ + 7)(2^ + 9) 
l)(£ -2)(^-3)(£-4)(£-5) 



(2£ - 11)(2£ - 7)(2^ - b){2i - 3){2£ - 1){2£ + 1){2£ + 3) 
21£{e -1){£- 2){£'^ - - 12^2 + + 54) 



3t-i{x) + 
3i-2{x) - 
3t{x) + 
je+2{x) - 
3e+A{x) + 
3e+6{x) 

je-7{x) - 
je-5{x) + 



{2£ - 9) {2£ - 7) (2^ - 3) (2^ - 1) (2^ + 1) {2£ + 3) {2£ + 5) 

35^(^2 _6)(-^4_;^4£2_^ 22) 
{2e - 1){2£ - 5)(2£ - 3)(2^ + 1)(2^ + 3)(2<' + b){2t + 7) 

35(£ + 1)(£2 + 2^ - 5) * + 4^3 _ g£2 „ 24£ + 9) 
(2^ - 5) (2f - 3) (2^ - 1 ) (2£ + 1) (2£ + 5) (2f + 7) (2f + 9) 
21 (£ +!){£ + 2){t + 3)(^'' + + q(2 _ 40^ + 15) 



(2f - 2,){2£ - l){2£ + 1)(2£ + 3)(2£ + 5)(2^ + 9)(2^ + 11) 



3i-s{x) - 
je-i{x) + 
ji+i{x) - 

je+3{x) + 
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7(^+l)(^ + 2)(^ + 3)(^ + 4)(^ + 5)(^2_^6^-3) . , , 



{2i - l){2t + l){2t + 3)(2^ + b){2l + 7)(2£ + 9)(2£ + 13) 
+ + 2)(^ + 3)(£ + 4)(£ + 5)(£ + 6)(£ + 7) 



J^+7(a;) 



(2£ + 1)(2£ + 3)(2£ + b){2£ + 7)(2£ + 9)(2£ + 11)(2£ + 13) 

l{t - - 2)(^ - 3)(^ - 4)(^ - 5)(£ - 6)(£ - 7) 

(2£ - 13)(2£ - 11)(2^ - 9)(2£ - 7)(2£ - 5)(2£ - 3)(2^ - 1)(2^ + 1 

2){l ^ 3)(£ - - 5)(2£^ - 101 - 19) 
(2£ - 13)(2£ - 9)(2^ - 7)(2^ - 5)(2£ - 3)(2£ - l){2(. + 1)(2^ + 3) 

14^(£ - 1)(£ - 2)(€ - 3)(2^4 - 12£3 - 22£2 + 120^ + 173) 
(2£ - 11)(2£ - 9)(2£ - 5)(2£ - 3)(2£ - 1)(2£ + 1)(2£ + 3)(2^ + 5) 
2U{1 - l){2f' - - 49^4 + 108^3 _^ 371^2 _ 436^ _ 873) 



(2£ - 9)(2£ - 7)(2£ - b){2l - 1)(2^ + 1)(2£ + 3)(2£ + 5)(2£ + 7) 



h-2{x) 



35(2^8 ^ _ 44^6 _ ;Lgo^5 _^ 286^4 + 848^^ _ qq^(2 _ ^q^q^ _^ g;^^) 

(2^ - 7)(2£ - 5)(2^ - 3)(2£ - 1)(2£ + 3)(2£ + 5)(2^ + 7)(2^ + 9) 
28(^ + 1)(^ + 2){2l^ + ld>t + 11^"* - 204^3 _ i^-jf ^ ggo£ _ 225) . 



(2^ - 5)(2^ - 3)(2£ - 1)(2£ + 1)(2£ + 3)(2^ + 7)(2^ + 9)(2£ + 11) 

IA{1 + + 2){i + 3)(£ + 4)(2£4 + 20^^ ^ 26^^ ~ 120^ + 45) 
(2£ - 3)(2£ - 1)(2£ + 1)(2^ + 3)(2£ + 5)(2£ + 7)(2£ + 11)(2£ + 13) 

4(^ + + 2)(^ + 3)(^ + 4)(£ + 5)(^ + 6)(2£^ + 14£ - 7) 
(2£ - 1)(2£ + 1)(2£ + 3)(2^ + 5)(2^ + 7)(2£ + 9)(2£ + 11)(2£ + 15) 
(£ + + 2)(^ + 3)(^ + 4)(^ + 5)(^ + 6)(^ + 7)(£ + 8) 



{2£ + 1)(2£ + 3)(2£ + 5)(2^ + 7)(2^ + 9)(2£ + 11)(2^ + 13)(2£ + 15) 
Alternative forms may by obtained by use of Eq. ( |All| ). 



jl+2{x) + 

■n+&{x) 4- 
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